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Abstract 

Using an SU{3) flavour symmetry breaking expansion in the quark 
mass, we determine the QCD component of the nucleon, Sigma and Xi 
mass splittings of the baryon octet due to up-down (and strange) quark 
mass differences in terms of the kaon mass splitting. Provided the average 
quark mass is kept constant, the expansion coefficients in our procedure can 
be determined from computationally cheaper simulations with mass degen- 
erate sea quarks and partially quenched valence quarks. Both the linear 
and quadratic terms in the SU{3) flavour symmetry breaking expansion 
are considered; it is found that the quadratic terms only change the result 
by a few percent, indicating that the expansion is highly convergent. 
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1 Introduction 



The masses of the baryon octet are now very accurately known, with results given 
in the Particle Data Group [T] as 

M^^P = 0.939565346(23) GeV, M^^^p = 0.938272013(23) GeV , 

M^"_P = 1.197449(30) GeV, M^^ = 1.18937(7) GeV , (1) 

M|^P = 1.32171(7) GeV, M^^ = 1.31486(20) GeV , 

around the outer ring of the octet and 

M^o^ = 1.192642(24) GeV, Mf^ = 1.115683(6) GeV, (2) 

at the centre. Isospin breaking effects (i.e. u - d quark mass differences and 
electromagnetic effects) are responsible for the nucleon, n - p, Sigma, S~ - S"*", 
and Xi, E~ - mass splittings 

(M„-Mp)""P = 1.293333(33) MeV, 
(Me--Ms+)""p = 8.079(76) MeV, (3) 
(Mh--Mho)^^p = 6.85(21) MeV. 

These are very small differences (particularly for the n-p mass splitting), ranging 
from about 0.15% to 0.7% of the baryon mass. (We shall not be considering the 
S° - A mass splitting here as these particles have the same quantum numbers and 
mix if isospin is violated.) In this article we shall be only looking at the hadronic 
or QCD contribution to these mass splittings, i.e. we are not going to consider 
electromagnetic effects. Both effects are small perturbations and can simply be 
added together. In the case of the n—p splitting we can argue that the hadronic 
effect is larger because the electromagnetic effect would tend to make the proton 
heavier than the neutron (as u quarks repel more than d quarks) which is not the 
case in the real world. There have been several previous lattice determinations 
of the QCD contribution to these mass splittings, e.g. [21E1I1], and also several 
lattice computations of the electromagnetic contribution, e.g. [21 [6l [TJ [U [9| [TO]. 
Non-lattice determinations include [llj . 

In Fig. [1] we sketch the lowest octet for the spin ^ baryons plotted in the I3-Y 
plane. The particles on the outer ring, namely the n{ddu), p{uud), T,^{dds), 
T,~{uus) and E~{ssd), H°(ssu) all consist of combinations of aab quarks (where 
we shall use the notation of denoting a quark, q, by a, 6, . . . which can be the up u, 
down d or strange s quark). In this notation a are the flavour doubly represented 
quarks, while b is the flavour singly represented quark. At the centre of the octet 
we have two states, A{uds) and S°(u(is), with the same quark content - u, d and 
s, but different wavef unctions. 

In [T3] we described a method for extrapolating from the SU{3) flavour 
symmetric point (where we have three mass-degenerate quarks) to the physical 
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Figure 1: The lowest octet for the spin 2 baryons plotted in the I3-Y plane. 

point, keeping the average of the quark masses constant. The form of the SU{3) 
flavour symmetry breaking expansion was developed both for non-degenerate u, d 
quark masses and for degenerate m, d quark masses. In [13] numerical simulations 
were performed for 2 + 1 flavours, i.e. with two degenerate light quark masses. 
Thus, effectively the mass 'average' of the ra, p baryon, and the S"*", S~ and 
S^, "EP baryons were considered. However, as the coefficients of the quark mass 
flavour symmetry are just functions of the average quark mass, the expansion 
coefficients do not change from using non-degenerate to 2 or 3 mass degenerate 
quarks provided that the average quark mass is kept constant p!3] . This gives us 
the opportunity to investigate isospin splittings, i.e. when the u quark mass is 
different to the d quark mass, using only results from 2+1 or 3 flavour simulations. 

As the baryon mass differences (e.g. n- p) depend on the u - d mass difference 
and are thus small, we flnd that it is sufficient to consider the SU{3) flavour 
symmetry breaking expansion in the quark mass including both linear terms 
(leading order, LO) and quadratic terms (next to leading order, NLO). These LO 
and NLO terms were given in |13j . 

We saw little curvature in hadron masses in the quark mass range considered 
in [13] (see also section [6]) , namely, from the three degenerate flavour pion mass 
at ~ 411 MeV to the physical pion mass ~ 140 MeV, so we conclude that a much 
larger quark mass range is needed to reliably determine curvature. We achieve 
this larger range by extending the numerical results to ^partially quenched^ or PQ 
quark masses (where the valence quark masses do not have to be the same as the 
sea or dynamical quark masses) with a spread of quark masses from about one 
third of the strange quark mass up to the charmed quark mass. We also consider 
part of the next to next to leading order (NNLO) (cubic terms) in the SU{3) 
flavour symmetry expansion. (The NNLO terms were also indicated in [TB]; we 
have now completed this computation, [2].) Thus we can consider a large quark 
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mass range to be able to determine the NLO or quadratic terms more accurately, 
using; the NNLO terms 'control'. 

Chiral perturbation theory looks at the breaking of the chiral SU{3) group 
by the quark masses - its expansion parameter is the quark mass itself, or equiv- 
alently, the masses of the pseudoscalar mesons. We are following an older strand, 
going back to Gell-Mann and Okubo [151 IB] of looking at the breaking of the 
non-chiral SU{3) symmetry, by quark mass differences. In chiral perturbation 
theory the expansion is about the point where all three quarks are massless; here 
we expand about a point where the three quarks have equal (non-zero) masses 
each about a third of the physical strange quark mass. 

As well as the baryon SU{3) octet flavour expansion, we will also need the 
values of the quark mass corresponding to the physical point. We can achieve 
this by considering the equivalent SU (3) flavour expansion, but now for the pseu- 
doscalar meson octet. The same procedure as for the baryon octet is required: 
first the SU{3) flavour expansion coefficients must be determined and then the 
experimental values of the masses of the K^, K'^ and mesons can be used 
to determine the required physical quark mass point. These can then be used, 
together with the 5*^7(3) baryon octet flavour expansion, to determine the mass 
splittings for the baryon octet. 

We shall find that the LO term is dominant (both for the baryon and pseu- 
doscalar octets) and so the NLO (and NNLO corrections) may be taken as an 
indication that our SU{3) flavour symmetry breaking expansion appears to be a 
highly convergent series. (This point is further discussed in Appendix lAl) 



2 Octet Baryons 

Before discussing partial quenching, we first consider the case where the valence 
quark masses are the same as the sea quark masses, the so-called ^unitary line\ 
The SU{3) flavour symmetry breaking expansion, [TB], for all of the outer ring 
octet baryons consisting of a pair of identical flavour quarks and a third, different 
quark can be compactly written up to NLO as 

M'^{aab) = + Ai{25ma + drrib) + A2{Smh - 5ma) 

+ Bol{Sml + 5ml + 5ml) (4) 
+ Bi{25ml + 5ml) + B2{5ml - 5ml) + B^{5mb - 5ma)'^ , 

with quarks q = a, b, . . . from (m, d, s), where 

5mq = mg — m, m = ^{m^ + m^ + mg) . (5) 

We shall consider the SU{3) flavour symmetry breaking expansion of M^{aab) 
\n\ . rather than M{aab). Of course from the viewpoint of the SU{3) flavour 
symmetry breaking expansion any function f{M) could be considered. For fitting 



4 



over a small quark mass range, a linear function is sufficient; for the large quark 
mass range considered in section 13.21 a better fit to the numerical data was found 
using M^(aa6) rather than M[aab). 

Note that Mq , Ai, A2, Bq, ...,-83 all depend on the average quark mass m, 
which will be held constant in the following simulations. Keeping m constant 
reduces the number of coefficients that must be determined (and indeed makes 
the computation tractable). 

From Fig. [1] we see that as there are six different masses on the baryon outer 
ring but just two linear parameters in LO, the fits are highly constrained. At 
the next order, NLO, in eq. (jl]) we are allowed four coefficients for the quadratic 
terms. 

We have in addition the trivial constraint 

6mu + Srud + 6ms = , (6) 

so we can eliminate one of these quantities if we wish to. 

Thus, to determine the octet baryon masses, we first have to determine the 
expansion coefficients and second we need to know the physical quark masses. 
In the following we shall denote the physical point by a star *. We thus have 
two distinct computations. As we shall see, the determination of the coefficients 
is helped by PQ simulations, while Sm* can be found by considering equivalent 
expansions for the pseudoscalar meson octet. 

A further problem is that the scale must be determined. As discussed in [13] , 
flavour blind (or gluonic) quantities are suitable. We denote these generically by 
X. One useful type of flavour blind quantity can be considered as the 'centre of 
mass' of the multiplet. Thus for the baryon octet, one possibility i^ 

= i(M2 + M^ + M|+ + Mi_ + M|o + M|_) . (7) 

At the physical point, from eq. ([T]), this gives 

X^"P = 1.1610 GeV. (8) 

In general for the SU{3) flavour breaking symmetry expansion we have from 
eq. ©, 

= Ml + 0{5ml). (9) 

Upon adding the masses the A2 and B2 terms vanish; while the Ai term vanishes 
upon using the constraint equation, eq ([6]), and thus this leads to the vanishing 
of the linear term in the quark mass. (This is indeed true for all flavour blind 
quantities.) 

^Another independent possibility would be X\ = ^(Af^ + M|n)- 
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Scale independent quantities can now be constructed by considering the ratio 
M"^ {aab) / Xjj . Expanding to NLO order in the quark mass, eq. (jl]) retains the 
same pattern and becomes 

(aa6) = I + Ai {26ma + 5mb) + A2 {Snib - Srria) 

- {Bi + 53) {5ml + 5rnl + 5ml) (10) 
+ Bii25ml + 5ml) + B2{5ml - 5ml) + B^{5mi, - 5m„)2 , 

where a ~ on a hadron mass means that it has been divided by X^r, so M = M/ X^v 
while a ~ on the expansion coefficients means that they have been divided by Mq , 
for example Ai = Ai/Mq. From eq. (Q we see that we have effectively replaced 
the |i?o term by —{Bi + B^). (However, this will not be important in the case 
discussed here; as for mass differences, these terms cancel again.) 
Alternatively we can re-write eq. ( ITOl) as 

M{aab) = 1 + A[{25ma + 5mb) + A2{5mb - 5ma) 

-liB, + B,){5ml + 5ml + 5ml) (11) 
+B[{25ml + 5ml) + S^(5m^ - 5ml) + B'^{5mb - 5maf , 

(essentially the equivalent SU{'i) flavour symmetry breaking expansion for M 



rather than M^) with 








= Ml, 




A'2 






B[ 


= UBi-lAl), 




B'2 


= i(52-|ili2 


), 


B', 


= i(53 + i(2ii 


-A2){A, + A2)). 


Although eq. ^ 


looks complicated. 


we shall only be 



(12) 



differences, which simplify the expressions. Writing the flavour expansions as a 
Taylor series in ± 5mu we find to NLO 



together with 
Me- - Ms+ 



and 



M{ddu) - M{uud) (13) 
{5md - 5mu) Ui - 2A'2 + {B[ - 25^) (^m^ + (5m„) 



M{dds) - M{uus) (14) 
{5md - 5mu) fei; - A'^ + {2B[ - 5^ + 2,B'^){5md + (5m„) 



M{ssd) - M{ssu) (15) 
{5md - 5mu) U'l + i'a + {B[ + B2 + 3&^){5md + 5m„) 
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where the constraint in eq. has been used to ehminate SitLs. These equations 
represent the different types of isospin differences possible for the baryon octet 
and are vahd over a range of quark masses. We shall consider them from the 
flavour symmetric point down to the physical point, determined by 5m*^—5m*^ and 
5m^ + (5m* = —5m*, 5m* being the physical point (which has to be determined). 

We have a check of these formulas, since the Coleman- Glashow relation [T8] 
should hold at this order, 

(M„ - Mp) - (Ms- - Ms+) + (Mh- - Mho) = 0(5mJ) , (16) 

which is indeed satisfied by eqs. f|T3|) - f|T5|l . 

We can find relations between isospin violation caused by m^ — m„ and to the 
SU{3) violation caused by m^ — mu, rris — rrid if we make any 5*3 transformation 
that changes the strange quark into a light quark and vice versa (for example 
■u— i-s— i-c/orrf-H-s); see Fig. [2l {Ss is the symmetry group of the equilateral 



< > 




Figure 2: Permutation group transformations link the isospin violation caused by 
md — rriu to the SU (3) violation caused by — m„, ms—rrid. The neutron-proton mass 
difference, Mn—Mp, dashed line, is mapped to M=o —My,+ ; the S splitting, Mj^- —Mj^+ , 
maps to M=- — Mp dot-dashed line; and M=- — M^o is related to M^- — M„, full line. 



triangle and a subgroup of SU{3).) Applying the transformation d ^ s on both 
sides of eqs. dH-dlSD we find to NLO 



My;- - Mn = {6ms-6mu)\A[ + A'2 + {B[ + B'2 + 3B'^){6ms + Sm^ 



(17) 



3 Determining the expansion coefficients 



We now first find the Ai, A^ and B\ 



B-s coefficients. 



7 



3.1 Partially quenched octet baryons 



Let us now generalise the previous results to the case when the valence quarks 
do not have to have the same mass as the sea quarks (i.e. we leave the unitary 
line). We have also generalised eq. (65) of [13] from NLO to NNLO [11]. We find 

(aab) = + Ai {2SfXa + 5fXb) + ^2 (SfXb - 5^a) 
+ lBo{6ml + 6ml + 6ml) 

+ Bi{26fil + 6fil) + B2{6fil - 6fj.l) + 3^(6 fit - 6iiaf 
+ Co6mu6md6ms 

+ [Ci{26fia + 6fib) + C2{6fib - Sfia)] {Sml + 6ml + 6ml) 

+ C3{6fla + 6flb)^ + Ci{6^a + 6flb)^{6fla - 6flb) 

+ C5{6fia + 6fib){6fia - Sfib)^ + Ce{6fla - 6flb)^ , (18) 

where a, b, . . . now denote three valence quarks of arbitrary mass, and we have 
defined 

^/^g = f^q-m qe{a,b,...}, (19) 

where fiq is the valence quark mass. In distinction to the sea quarks, there is 
no restriction of the form of eq. ([6]) on the values of valence quark masses. The 
numerical values of the Mq, Ai, A2, Bq, Bi, . . ., B^ and Cq, . . ., Cq coefficients 
are the same for PQ as for the unitary case. 

While we see that this is a relatively straightforward generalisation of eq. (jl]), 
we note that the term proportional to Bq remains unchanged. In addition the 
Cq term also depends entirely on sea terms, while the Ci and C2 terms are a 
mixture of sea and valence terms. Thus if we wish to determine these coefficients 
we must vary the sea quark masses; to determine the other coefficients it is 
sufficient to vary the valence quark masses alone, while keeping the sea quark 
masses constant. So this gives the possibility of extending the (computationally 
expensive) sea quark mass simulations with (computationally cheaper) valence 
quark mass simulations to determine most of the coefficients. If we work on a 
single sea background, then the Cq term can be absorbed into the Mq term, while 
the Ci and C2 terms can be absorbed into the Ai and A2 terms. If we vary the 
sea quark masses this allows a determination of these coefficients. However due 
to the constraint m = const., or equivalently eq. (E]), 6mq cannot vary much and 
we know from [13] that in this range the LO dominates, so these coefficients are 
difficult to determine and contribute just noise. So practically we shall ignore 
these terms in fits. (Alternatively, the constraint m = const, could be relaxed, 
but then we have additional expansion coefficients, which we wish to avoid.) Thus 
in this article we regard the NNLO terms as 'control' on the LO and NLO terms. 

As discussed in section [2l we can consider scale independent quantities. Thus 
in analogy to eq. f lTTj) we have 

M^{aab) = 1 + Ai{26fXa + 6fib) + M^f^b - 6fXa) 
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- {Bi + B^){5ml + 5ml + 5ml) 

+ 5i(25/i2 + 5^il) + B2{5^il - 5fil) + Bs{5fi, - 5^iaf 

+ ((73 - ?>Cr^)5mu5md5ms 

+ Ci{25^a + 5^b) + C2{5^JLb - 5^a) i^'rnl + S'TiT'd + ^^l) 
+ C'5(5/ia + 5nb){5iJ^a - 5fXbf + C6i5fia - 5^lbf , 



(20) 



(where Xn always depends just on the sea quarks and is given by the NNLO 
extension of eq. Q). 

Furthermore, these equations remain vahd if two of the sea quarks are degen- 
erate in mass, i.e. m„ = m^ = mi, the crucial point being that m must remain 
constant (as all the coefficients are functions of fn). This means that from dynam- 
ical 2-1-1 flavour simulations we can determine the u-d mass splittings. The only 
change to eq. (!20!) when m„ = m^ is that some terms become slightly simpler. 



5ml + + ~^ 65ml , 5mu5ma5ms — )■ —25mf , 



where we have used the constraint equation, ([ 

5ms = —25mi . 



which now becomes 



(21) 
(22) 



This gives 

M^(aa6) 



1 + Ai{25fla + 5 fib) + A2{5flb - 5fla) 

- 6{Bi + B^)5ml 

-f Si(25/i2 + 5^il) + B2{5^il - 5fil) + Bs{5fib - 5fi,)^ 

- 2(^3 - 3C5)5mf 



+ 6 



Ci{25fla + 5fib) + C2{5flb - 5fia) 



5ml 



+ C3{5fia + 5fib)^ + C4,{5fia + 5fib)^{5fla " 5Hb 
+ C5{5fla + 5flb){5fia - 5flb)^ + CG{5fia " 5 fib) 



with 



UM^IU) + M\Us) + M^ssl)) . 



(23) 



(24) 



(For a quark mass-degenerate 3 flavour simulation eq. 0231) simplifies further as 
5mi = = 5ms.) In other words, using eq. (125]) gives us all the information we 
need to find the quark mass contribution relevant for the 1 -|- 1 -|- 1 case. 



3.2 Numerical results 

Simulations have been performed using Nf = 2-1-1 0(a) improved clover fermions 
[19] at /3 = 5.50 and on 32^ x 64 lattice sizes, as described in more detail in 
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j. Errors given here are statistical (using ~ 0(1500) configurations); possible 
systematic errors are discussed in Appendix |X] and incorporated into the final 
results in section [3 

A particular starting value for the degenerate sea quark mass, mo, is chosen 
on the SU{3) fiavour symmetric line, and the subsequent sea quark mass points 
mi, rris have then been arranged in the various simulations to have constant fn 
(= mo). This ensures that the expansion coefficients do not change. It was found 
in [13] that a linear fit provides a good description of the numerical data over 
the relatively short distance from the symmetric point down to the physical pion 
mass. This helped us in choosing the initial point on the SU{3) fiavour symmetric 
line to give a path that hits (or is very close to) the physical point. 

In a little more detail, the bare quark masses in lattice units are defined as 

mg = l-(- —] with qe{l,s,0}, (25) 



2 \Kq Kq-c 



(together with eq. (|5]) for 6mg) with the index g = denoting the common mass 
degenerate quarks along the SU (3) fiavour symmetric line, and where vanishing 
of the quark mass along this line determines ko-c- Keeping m = constant = mo 
gives 

(26) 



_3 _2_ 



So once we decide on a ki this then determines k^. Note that kq-c drops out of 
eq. ( 126|) . so we do not need its explicit value. The initial SU{3) fiavour symmetric 
kq value chosen here, namely kq = 0.12090, [13] is very close to a point on the 
path that leads to the physical point. The constancy of fiavour singlet quantities 
over the range from the SU{3) fiavour symmetric line down to the physical point 
|13] . leads directly from X^^ to an estimate for the pion mass of ~ 411 MeV (i.e. 
eq. ( 14 7p ) and similarly from X^v a value of the lattice spacing of a ~ 0.079 fm. 

While, as discussed earlier, simulations between the SU (3) fiavour symmetric 
point and the physical point are in principle enough to determine the linear and 
quadratic expansion coefficients, in practice the range is not sufficiently large 
to reliably determine the quadratic terms. In order to determine the quadratic 
coefficients more precisely, additional PQ simulations have been performed on the 
set of gauge configurations that have all three sea quark masses equal, i.e. at the 
SU{3) fiavour symmetric point kq = 0.12090. For these particular simulations 
6mi = = 6ms automatically, fig is defined identically to m^, eq ( 125|) . by 
replacing m^ — )• fig with g G {a,b, . . . ,} together with eq. ( IT^ for 6 fig so that 

Sfig = fig-m, fi^ = l-(- —\. (27) 

I \Kg KQ-C J 

We have chosen a wide range of PQ masses starting from a slightly heavier mass 
than mo (to avoid any possibility of so-called 'exceptional configurations') and 
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reaching up to masses ~ rricharm- The values are given in Appendix [B] in Table HI 
A direct fit is made to these PQ masses, the unitary masses of pj^j, which all 
have the same fixed m, and additionally three PQ masses, which we call Mjy^ = 
M{sss) again with the same fixed m. (This is analogous to the pseudoscalar t]s 
considered later.) These values are also given in Appendix [B] in Tabled Using 
the fit function of eq. fl20|) and ignoring the Ci and C2 terms as discussed in 
section [3?T] gives the results in Tabled! together with the NNLO coefficient values 





A2 


Bi 


B2 


B3 


10.15(12) 


1.828(157) 


13.51(126) 


-10.29(139) 


-14.90(144) 



Table 1: Results for the baryon octet expansion coefficients. 

of C3 = -4.60(115), = -17.5(26), C5 = -1.88(310) and Cq = -3.65(184). 
We note that in comparison to the NLO coefficients, the NNLO coefficients are 
poorly determined. The bootstrap (MINUIT) fit used gave x^/dof ~ 0.4. 

We now compare these results to a plot for illustration. The simplest to 
consider is setting 6fia = 5/^6, i.e. degenerate valence quark masses. For simplicity, 
but slightly inaccurately, we shall in the following simply say b = a. (Of course 
we still need two different quark flavours.) With 6mi = 0, eq. (123|) then becomes 

= A, + B,6fia + IC36fil. 28 

In Fig. [3] we plot {M'^{aaa) — l)/(35/ia) against 6fia for the PQ data (together 
with the cubic fit coefficients from Table [T]) which due to the denominator is a 
sensitive plot. There is good agreement. (We postpone the comparison to the 
unitary data, 'fan' plots, until section [61 ) 

4 Octet pseudoscalar mesons 

Determining the octet pseudoscalar mass splittings (or more accurately the split- 
tings of the quadratic masses) will give 5m*, 6m^ and 5m* (the quark masses 
at the physical point). This closely follows the baryon octet procedure, we must 
again consider the analogous flavour symmetry expansion for the pseudoscalar 
meson octet together with the known experimental masses of the pions and kaons. 

In Fig. m we sketch the lowest pseudoscalar octet in the I3-Y plane. We 
have K'^{ds), K~^{us), n^{ud) together with K (sd), K~{su) and 7i~{du) in the 
(outer) ring of the octet. From charge conjugation or C invariance we further 
have Mj^o = M^o, M^- = Mk+, and M^- = Mn-+ (which is in distinction to the 
baryon octet which does not have this constraint). 
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4.1 PQ pseudoscalar meson flavour expansions 

The corresponding formulas for the octet pseudoscalar mesons are simpler than 
for the octet baryon, due to the constraints imposed by C invariance. The follow- 
ing S'f/(3) flavour breaking expansion formula is always valid for quarks q = a, 
6, . . . , in (m, rf, s, . . . , ) 

M^{ab) = M^^ + a{Sma + 5mb) (29) 
+ (3ol{6ml + 6ml + 6ml) + /3i{6ml + 6ml) + /32{6ma - 6mb)^ , 

in quark masses up to the NLO as discussed in [13j. (Note that M{ab) = M{ha).) 

Again combinations of masses can be chosen, so that due to eq. ([6]) the linear 
term in eq. ( 129|) vanishes, which is equivalent to averaging the outer ring of 
particles and finding the 'centre of mass' of the octet. In particular if we set 

Xl = l[Ml+ + Mlo + Ml+ + Ml- + + Ml-) 

= \{Ml^ + Ml, + Ml^) , (30) 

this gives 

Xl = Ml+{lPo + lPi + P2)i6ml + 6ml + 6ml) 

= Ml + 0{6ml). (31) 

In the partially quenched case, the above SU{3) flavour expansion can be 
generalised to 

M^{ab) = M^^ + a{6fia + 6fib) 
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Figure 4: The lowest octet for the spin pseudoscalar mesons plotted in the I^-Y 
plane. 

+ 'yo5mu6md6ms + 7i (5/ia + SfJ'b)i6ml + 5ml + 6m1) 

+ 72('5/ia + SjJbf + Isi^l^a + 5lJb){5lJ^a " S fXhY , (32) 

where the NLO was also discussed in [13] and again we have extended the formula 
to the NNLO case [l^ . This is again a general formula valid for possibly differing 
masses of the sea, m^, and valence quarks, /iq. The same notation has been 
used as in sections [2] and El (In particular remember that is unconstrained 
in distinction to the sea quarks, eq. ([6]). The unitary line is recovered when 
liq mg.) Again as with the PQ octet baryon case, eq. ( !T8|) . we see that at the 
NNLO there is a term, the 71 term, which mixes the sea and valence quarks; as 
discussed in section 13.11 this again makes the numerical determination of these 
coefficients difficult (we cannot vary the sea quark masses over a large enough 
range). So in the same spirit as section 1371] we ignore this term and regard the 
NNLO terms as 'control' terms. 

If we wish to consider 'physical ratios' then the masses can again be conve- 
niently normalised using X^. Expanding to NNLO, we find 

M\ab) = 1 + a{5fXa + 5f^b) 

- (lA + ^2){Sml + 6ml + 6ml) + Pii^l^l + ^f^l) + P2{6fia - 6^ibf 
+ 2(72 - 3%)6mu6md6ms + ^i{6fia + 6^,b){6ml + 6ml + 6ml) 

+ 72(5yUa + 6^.bf + 73(5yUa + 6flb){6fla - 6^bf , (33) 

where again a ~ on a hadron mass squared means that it has been divided by 
Xl (which only depends on the sea quarks) while on an expansion coefficient it 
means that the coefficient has been divided by Mq^ for example a = a/M^^. 
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Again for a 2 + 1 flavour simulation (the case that will be considered here) 
eq. ( 133|) slightly simplifies to become 

M'^{ab) = 1 + a{5fXa + 5 fib) 

- 2(2/3i + 3/32)5m2 + ^,{6fil + 6fil) + hi^fi^ - 6fn,f 

- 4(72 - 373)(5mf + 671 (5/^^ + 5fXb)Smf 

+ l2{Sfla + Sfibf + ^siSfla + 6fib){Sfia " Sflbf , (34) 

with the same coefficients, provided of course that m remains constant and where 

X'^ = l{2M\ls) + M^fl)) . (35) 

4.2 Determination of the coefficients 

As in section |3] we need to consider PQ masses at the SU (3) flavour symmetric 
point for the sea quark masses (see Table jointly with the 32^ x 64 lattice size 
unitary results from Table 20 of [13] all with the same rn constant value. From 
eq. we need to determine the constants a, (3i, (32 and 72, 73. 

A 5-parameter fit to the PQ and unitary data then yields the results of Table [2] 
together with the NNLO coefficient values of 72 = -16.4(12) and 73 = -20.3(39). 



a 


Pi 


^2 


41.17(8) 


76.50(148) 


-45.81(189) 



Table 2: Results for the pseudoscalar octet expansion coefficients. 

The bootstrap (MINUIT) fit used gave xV^of ~ 1.7. 

As in section [3^ it is useful to compare these fit results in a plot. In parallel 
to eq. (1251) let us again consider the simple case of degenerate quark mass, i.e. 
where quark and antiquark have the same mass, but different flavours, so they 
cannot annihilate. We set 6 = a and Smi = in eq. (IMl) and so consider 

= a + + 4725/x^ (36) 

In Fig. |5]we plot (M^(aa) — l)/(25/ia) against 6fia using the PQ data. This is 
compared with the cubic fit from eq. f l36|) and Table [2j There is good agreement. 

4.3 Pseudoscalar meson isospin spfittings 

Having determined a, Pi, P2 we can now find 5m„, 6md, Srris given the masses 
around the outer ring of the pseudoscalar octet. We postpone a discussion of 
the numerical values until the next section and here just derive the relevant 
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Figure 5: (M^(aa) — l)/(26fia) against 6fia- Tlie full line is eq. ([36|) . using the a, f3i 
fit values from Tableland 72 (given in the text). 

formulas. As they are trivially valid over a range of quark masses (and not just 
at the physical point) we give these more general expressions here. 

From eqs. (fT3l)-(|T5l) we see that we need Srud — Srriu = and Snid + 5mu = 
6m . So as in section [2] we consider the mass difference 

MIo-Mk+ = {6md - 6m^)[& + 01 + 3(32) {Smd + Sm^)] 

= Sm^a + 01 + 3(32)5m+] . (37) 

Alternatively we can consider 

M^K+-M^+ = i^ms - 6md)[a + 01 + 3f32){6ms + 6md)] 

= -i(5m_+35m+)[5 + i(/3i + 3/32)(5m_-5m+)], (38) 

and 

Ml,-Ml+ = {5ms-5mu)[a + 0i + 3(32)i5ms + 5mu)] 

= i(5m_-35m+)[a-i(/3i + 3/32)(5m_ + 5m+)], (39) 

where we have used the constraint, eq. (E]), to eliminate Srus in the second line 
of eqs. ( 138|) and ( l39l) to re- write the equations in terms of 5m_, (5m+. 

Of course only two of the equations ( 137|) -( |39|) are independent. Choosing 
eqs. (137|) and ( |38|) . these quadratic equations can be solved iteratively to give 
SrUzf. We start the iteration from the linear term alone or LO in the quark mass. 
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From eq. ( 137|) we see that is sufficient to determine 5m+ to LO (and (5m_ to 
NLO). Thus we find 



2 



Sma + Sruu = -— (K^^ko + M^k+) ' M^+) [1 + ^nlo] • (40) 

As 2(/3i + 3/32)/3a;^ ~ —0.02 then as, in particular, at the physical point ^(M^o + 
M^%) - M*+ is ~ 1, we expect that the NLO term for Srri- is smalljl i.e ~ 3% 
and tends to reduce the value of 6m- slightly. (See next section for the numerical 
values.) 

5 Physical values of the quark masses 



To proceed further we now need to substitute eqs. fHOj) into eqs. f[T3l) - f[T5|) to 
give the pure QCD contribution to baryon octet mass splittings in terms of the 
pseudoscalar octet masses. 

Before considering this however (see section EH), we shall first discuss electro- 
magnetic effects and determine the physical values of the quark masses (5m* and 
(5m^ given the experimental values of the pseudoscalar masses. This will enable 
us to investigate the convergence of the SU{3) fiavour breaking expansion. The 
experimental masses are, pL], 

M^+P = 0.13957018(35) GeV, 
M^j^^ = 0.497672(31) GeV, 

M^"P = 0.493677(16) GeV, (41) 

(with as already mentioned M-^ = M^o, M^- = Mk+ and M^- = M^+), on the 
outer octet ring, and at the centre 

M^"P = 0.1349766(6) GeV. (42) 

We now need to consider electromagnetic effects (which may now be compa- 
rable to the u - d quark mass difference which is also small). Electromagnetic 
effects tend to increase the mass of charged particles (due to the photon cloud). 



^This is also true for 5m+. The NLO term is 



3(5 



,(M|o+Af^+)-M2, 



Again together with (/3i + 3l32)/3a^ 0.01 and ^{M'^l + M*^\) - M*^ 1, this means that 

the correction NLO term for 5ra+ is also small. 
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As a help to estimate these unknown effects, we use Dashen's theorem, [20], which 
states that if electromagnetic effects are the only source of breaking of isospin 
symmetry (i.e. m„ = m^), the leading electromagnetic energy contribution to the 
neutral pseudoscalar particles, i.e. the 7r°, K^, vanishes, while that due to the 
charged particles, i.e. the vr^, is equal. As the mass difference in and 
due to the u - d quark mass difference is neghgible, 0(0.1 MeV), e.g. [17], we can 
thus write (e.g. [^[22]) 

M^^r = M*^l+^,,, M'^^' = M*^l. (43) 

(The * denotes values at the physical point for the pure QCD case.) Dashen's the- 
orem has corrections of 0{aQEDTTtq) from higher order terms. Sometimes possible 
violations of the theorem are parametrised by [H 

M^^o - M*^l = (Mio - M^,)^^"^ + (1 + e,) (M^^ - M^^o)''^^ , (44) 

where = corresponds to Dashen's theorem. For example a positive value for 
would thus tend to increase slightly the value of M^l — M^+. From the first 
equation in eq. (HOjl this would only affect the leading term in our expansion, so 
our main result in section 17.21 will be given in terms of the kaon mass splitting, 
M^o — and where we shall consider as an additional error. 
Inverting the relations in eq. (l43ll gives [1] 

M*J = M17\ M;,l=M'^l'-iM:r'-Ml7'), M*^l = M'^r, (45) 

or 

M*+ = 0. 13498 GeV , = 0.49240 GeV , M^o = 0.49767 GeV , (46) 

which we shall use as the pure QCD values due to differences in the quark masses 
with the electromagnetic effects subtracted away (assuming Dashen's theorem). 
This gives from eq. (IHUIl l^. 

= 0.4116 GeV, (47) 

(of course this is very close to the experimental value of X^^^ = 0.4126 GeV). 

In Table[3]using the masses given in eq. fH6l) we give our results, first giving the 
LO results for the quark masses, then the next line, NLO", gives the results from 
eq. fHOl) . which we will be using here. As a check, the third line compares this NLO 
result to the NLO result using the Newton-Raphson method to iterate eqs. (1571) 
and (155]) . We find the Newton-Raphson procedure converges very quickly (after 

•^In [13] we used the average kaon and pion masses as we were strictly in the 2 + 1 flavour 
case. Here we need to take into account the differences between charged and neutral mesons. 
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5m*^ — Sml 


Sm*^ + Sml 


Sml 


Sm*^ 


Sm* 


LO 

NLC^ 


0.0007485(14) 
0.0007245(27) 
0.0007249 


-0.02168(4) 
-0.02204(6) 
-0.02204 


-0.01121(2) 
-0.01138(3) 
-0.01138 


-0.01047(2) 
-0.01066(3) 
-0.01066 


0.02168(4) 
0.02204(6) 
0.02204 



Table 3: Results for the bare quark mass in lattice units at the physical point using 
eq. dMI) as input. The first line is the LO result. The next line, NLO*^, uses the result 
of eq. (|40p . while as a check the last line in italics, NLO^, iterates eqs. (|37p and (|38p . 

one or two iterations) and find good agreement between the two results (well 
within the error bars of NLO" in Table |3]) and so we can be confident that 
eq. ( HOj) is a good approximation for the inversion of eqs. (137|) and (138|) . 

We also note that the differences between the LO and NLO are small of the 
order of a few percent, indicating that the SU{3) flavour expansion about the 
flavour symmetric point appears to be highly convergent. (We shall discuss this 
a little further in section |6]).) 

6 Comparison with 'fan' plots 

We now compare the fit results with the mass values from previous 'fan' plots, 
[13] . which describe the evolution of the pseudoscalar and baryon octet masses 
along a path from the SU{3) symmetric point down to the physical point. In 
[13] the isospin degenerate limit, i.e. = ^ mi, was considered. So for this 
comparison we take the physical quark mass, in lattice units, as 

Sm; = ^{Sml + Sm^) = -0.01102(3) . (48) 

In Fig. [6] we show numerical results between the SU (3) flavour symmetric 
point and the 'physical point' for the numerical pseudoscalar octet on the unitary 
line (keeping m = const.) from [13]. These masses are compared to the quadratic 
flt using eq. (i.e. together with just the results of Table |2]) for the 2-1-1 flavour 
case, i.e. m„ = m^ = mi. The NNLO terms in the SU{3) flavour symmetric 
expansion can be safely ignored in the small Smi range. Compare the scale of the 
X-axis of Fig. [3] with that of Fig. [6l We consider M^(/I), Mk{Is) and the flctitious 
PQ particle rjs, with mass M(ss). The comparison is satisfactory. We also show 
results from [13j on smaller 24^ x 48 sized lattices. This shows that flnite size 
effects for these ratios are rather small. Of course, the value of Sm'^ just serves 
to deflne pseudoscalar meson mass ratios at the 'physical point' in a Nf = 2 -|- 1 
flavour world. For completeness we give these numbers here: M*^ = 0.1077(26), 
= 1.446(1) and M*^ = 2.884(4). 

In Fig. [7] we show the comparable baryon octet results. As well as considering 
the nucleon mass, Mjv (///), we also show Mj]{lls), M-^{lls) and also a flctitious PQ 
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Figure 6: The pseudoscalar meson octet 'fan' plot, M'^^/X'^ {-kq = vr, K, rjs) versus 
6mi. The numerical mass values are taken from [13] where filled up triangles, left 
triangles, down triangles are '/r(//), K{ls), r]s{ss) results respectively using 32^ X 64 sized 
lattices. The common symmetric point is the filled circle. The opaque up triangles, 
left triangles, down triangles are from 24^ x 48 sized lattices (and not used in the fits 
here) . The quadratic fits are taken from eq. (j34p , together with Table [2j The vertical 
line from eq. (j48p is the pure Nj = 2 + 1 QCD physical point, while the opaque circles 
are the pure QCD hadron mass ratios for 2 + 1 quark flavours. 



particle - Mn^{sss). Again the comparison of the NLO (quadratic) fit, using the 
expansion coefficient values given in Table [H to the numerical data is satisfactory. 
For completeness we give here the values at the 2 + 1 QCD physical point (opaque 
circles) of M^^ ^ 0.6704(46), M^^ ^ i.051(6), Mi^ ^ 1.278(8) and M*/^ = 
1.692(9). For a comparison to these values, the stars in Fig. [7] represent the 
average of the squared masses of the appropriate particle, as defined in the figure 
caption. 



7 Results and discussion 

We are now in a position to numerically determine the baryon octet mass split- 
tings due to pure QCD u - d quark mass differences (see section I7.H our main 
result), and then in section W72\ estimate physical values for the splittings. 
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V: igure 7: The baryon octet 'fan' plot, Ml,^/X% {Nq = N, S, E, Ns) versus 6mi. 
The numerical mass values are taken from [13j where filled up-triangles, left-triangles, 
right-triangles, down-triangles are N{lll), T,{lls), E(ssl), Ns{sss) results respectively 
using 32^ X 64 sized lattices. The common symmetric point is the filled circle. The 
opaque up-triangles, left-triangles, right-triangles, down-triangles are from 24^ x 48 
sized lattices (and not used in the fits here). The quadratic fits are from eq. (j23p . 
together with Table [H The vertical line from eq. (j48|) is the Nj = 2 + 1 pure QCD 
physical point, with the opaque circles being the determined pure QCD hadron mass 
ratios for 2 + 1 quark flavours. For comparison, the stars represent the average of 
the squared masses of the appropriate particle on the outer ring of the baryon octet. 
Fig. m i.e. M;;^{m) = (Mr''^(cid«) + Mp""P^(TO(i))/2, M^'^{lls) = {M^^^{dds) + 
M^^ '^iuus))/2, Mi^ssl) = (M|"P \ssd) + M|^p \ssu))/2. 



7.1 Baryon octet mass splittings 

After re-writing quark masses in terms of the pseudoscalar octet masses, sec- 
tion 14.3^ and finding that the expansion is highly convergent in the relevant 
quark mass range, e.g. Table [3l we now insert the expansion of eq. (HO!) into 
eqs. (fT3l) - (fT5|) which gives the reasonably compact results 



Mat - Mjv/ = An_n' 
for the pairs {N, N') 



1 + Bn_n> [UMko + Ml^)- m: 
= {n,p), and where 



M; 



M' 



,(49) 



A. 



A[ — 2A2 



n—p 



B, 



a 



n—p 




h + 3/32 B[ - 2B', 



(50) 
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together with 



and 



2A[ — A2 



Bs 



a 




j3i + 3^2 -B'2 + 35^ 



2^1 — ^2 



(51) 




/3i + 3;g2 5; + 5^ + W'^ 



a 



A[+A'2 



(52) 



As discussed before, as we have unknown electromagnetic effects, then we 
shall first present our results in a general form with (known) coefficients between 
the baryon and pseudoscalar meson splittings, within a pure QCD context. We 
now insert the numerically determined values from Tables [1] and [2] into eq. fH9|) , 
together with eqs. (I50|) - (l52!) . to give 



M„ - M„ 



together with 



0.0789(41) (8) (8) (32) M|o - Mk+ 



1 + 0.0817(92) i\{Mlo + Ml+) - Ml 



(53) 



M^--M^+ = 0.2243(35) (22) (2) (90) (m|,o -M|-+ 



X 



1 + 0.0077(30) ( liMlo + Mk+) ' ^1+ 



(54) 



and 



0.1455(24)(13)(6)(58) [mI, - Ml+ 



X 



1 - 0.0324(50) 



ml. 



K+, 



Mi 



(55) 



as our main numerical result. This is a pure QCD result: the isospin breaking 
is just due to the different masses of the m, d and s quarks. (To recapitulate, * 
means at the physical point and for the baryon octet M„ = M^/Xtv, etc. and for 
the pseudoscalar octet M^o = M^o/Xjr, etc. where X'^ is the 'average' hadron 
(mass)^ of the 'outer' ring of the octet, given numerically in eqs. ([H]) and ( H71) .) 
Even with the long lever arm of PQ some of the NLO terms, i.e. for — S+ are 
only marginally determined. 

The first error is statistical; the other errors are systematic as discussed in 
Appendix |Al The second error bar is due to possible finite size effects, the third 
error estimates a possible error due to convergence of the SU (3) flavour breaking 
expansion, while the last error is due to the choice of path to the physical point. 

The above results are valid for a range of quark masses; we shall now specialise 
to the physical point, as discussed in section [51 
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7.2 Physical values of the mass splittings 

We first note that the NLO term is always small of the order of a few percent, 
and only slowly decreases for increasing number of strange quarks in the baryon. 
Using for the NLO term the mass values given in eq. (1461) for the pure QCD case 
we have ^(M^o + M^+) — M*+ = 1.339 (using the experimental rather than pure 
QCD masses would introduce negligible additional errors into the NLO term). In 
addition, using Xj^^ and X*^ from eqs. ([8]) and P7|) . this gives 

M*-M; = 0.599(32)(34)(M^^o-M^i), 
M*_-M*+ = L553(25)(63)(M^2^-M^i), 
Mi- - M5o = 0.954(18) (41) (M^^o - M*j/+) , (56) 

where all masses in these equations are now in GeV, as our final result in terms 
of the pseudoscalar kaon masses. The first error is statistical, while the second is 
the total systematic error. 

Using again the kaon values from eq. f H^ and regarding possible violations 
of Dashen's theorem, eq. (l44l) . as a further systematic error, our isospin breaking 
effects due to pure QCD alone (in MeV) are 

M*-M; = 3.13(15)(16)(76|eT,|)MeV, 
M*_-M*+ = 8.10(14)(33)(193|e^|)MeV, 
M|_-M|o = 4.98(10)(21)(120|e^|)MeV. (57) 

In general, comparing eq. fl57|) with eq. ([3]) would indicate that electromagnetic 
effects for n — p are negative, for — S+ are small and for — S° are positive. 

The uncertainty due to is the dominant error once we convert to MeV 
(i.e. unknown EM effects are the largest source of error). As an example, tak- 
ing = 0.7, H |23] gives M* - M; = 3.13(15)(16)(53) MeV, M*_ - = 
8.10(14)(33)(135) MeV and - M^o = 4.98(10)(21)(84) MeV. 

In |llj , a determination of the n — p isospin breaking effects due to electro- 
magnetic effects was given, with a result of —1.30(47) MeV. Adding this to the 
result of eq. gives 

(M„ - Mp)*+«^° = 1.83(52)(76|e^|)MeV. (58) 

This is to be compared with the experimental result given in eq. Q of 1.29 MeV. 
Thus we find consistency (within errors even with |e^| ~ 0). Thus this result also 
indicates that violations of Dashen's theorem seem to be small. 

In Fig. [8] we compare this n - p mass difference including QED, (M„ — 
Mp)*+'3^^, eq. dSHl) with = 0.7, together the results of |2l |3l i]. The filled 
circles use the QED determination of ^Jj, while the filled square includes the 
full determination from that reference. Despite the fact that QED effects are 
treated slightly differently in each work, good agreement amongst the various 
determinations and with the experimental result is found. 
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Figure 8: Comparison of the n - p mass difference of the present result (QCDSF- 
UKQCD or lowest plotted number) with NPLQCD, Blum et al., and RM123 [SlElll], 
respectively (top to bottom plotted numbers). The filled circles use the QED determi- 
nation of while the filled square gives the full QCD and QED determination from 
Blum et al. . The vertical dashed line is the experimental result given in eq. ([3]). 

7.3 Conclusions 

In this article we have introduced a method to determine the isospin breaking 
effects due to QCD for octet baryons. Using an SU (3) flavour symmetry breaking 
expansion in the quark masses, the pseudoscalar meson octet expansion coeffi- 
cients can be determined, which leads to an estimate at the physical quark mass 
point. This can then be used together with the equivalent SU{3) flavour sym- 
metry breaking expansion for the baryon octet to determine the baryon mass 
splittings. The expansion coefficients depend only on the average quark mass m. 
Thus we can use degenerate sea quark masses (either with mass degenerate u and 
d quarks or with mass degenerate u, d and s quarks, i.e. at the flavour symmetric 
point) provided that m remains unchanged. These are computationally cheaper 
simulations than those with mass non-degenerate m, d and s quarks. 

A further advantage of our procedure is that the mass expansion formulas 
are easily generahsed to the case of differing valence quark masses to sea quark 
masses (i.e. partially quenched valence quarks). This allows an extension of the 
quark mass range to heavier quark masses, so that both LO and NLO coefficients 
(i.e. linear terms and quadratic terms in quark masses) can be determined. Our 
final results are given in eq. flFTl) . As the NLO turns out to be only a small 
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Blum et al 



RM123 



QCDSF-UKQCD 
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correction to the LO this gives us confidence that the SU{3) fiavour symmetry 
expansion appears to be a highly convergent series. 
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Appendix 



A Systematic errors 

We now consider in this Appendix possible sources of systematic errors: finite 
lattice volume, convergence of the SU{3) flavour breaking expansion, the path to 
the physical point and finite lattice spacing. 

Finite lattice volume 

Comparing the available 24^ x 48 with the 32^ x 64 lattice data in Fig. Vindicates 
that for these mass ratios finite size effects are small. (Finite volume effects 
were also discussed in [13] in section 8.3.1.) We now use these unitary results 
to estimate possible finite size effects. For the lightest 24^ x 48 point, namely, 
= (0.121040, 0.120620),^M^L ~ 3.4 (where L is the spatial lattice size, 
here 24a), and the mass ratio M has finite size error ~ 1%, from comparing the 
24^ X 48 lattice result with the 32^ x 64 lattice result. (Either using the results of 
[13] directly or equivalently taking the square root of the results in Fig. [71) On the 
lightest 32^ x 64 lattice point, namely (0.121145,0.120413) we have M^L ~ 3.1, 
so we expect the finite size errors in the mass ratio Mn will be approximately 
the same, i.e. also ~ 1%. We use this in section [7] to estimate systematic errors 
arising from finite volume effects. 

SU{3) flavour breaking expansion 

We first note that in Fig. [7]in the range \Smi\ < 0.01 (and \6ms\ < 0.02), there is 
little curvature. This is in agreement with the SU{3) flavour breaking expansion, 
eq. (ITU]) , where each order is multiplied by a further Srriq. From Tabled] for the 
Ai and Bi coefficients (and for the NNLO order the Ci coefficients) we see there 
that they remain approximately all of the same order, so we expect that every 
increase in the order leads to a decrease by about an order of magnitude in the 
series. This is confirmed in the present case as we have compared the NNLO 
determination of Mat, Ms and M= with the NLO results (removing the heavier 
quark mass points until the x^/dof is approximately the same). The change in 
M^v was less than half a percent and for M^, M= less, which is equivalent to 
using ~ 10% of the NLO term to estimate systematic errors. 

A further example to illustrate this convergence is given by eq. (42) of [T3] . 
Here we can form sums and differences of the decuplet masses which are of order 
(5m°, 6ml, ^^"^ Srnf. We see that each time we add a factor of 5mi the 
quantity decreases by an order of magnitude (in fact usually by a factor of ~ 20), 
and the 0{Smf) quantity is about 2000 times smaller than the leading order 
quantity. So we believe that convergence is very good for hyperons. Such an 
expansion is very good compared to most approaches available to QCD. 
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Path to physical point 

As also discussed in [13], the chosen trajectory in the - mi plane (keeping 
m constant), appears to not quite go through the physical point. Using (see 
eq. (!35l) ) and Xn then from Table 15 of [13j, we see that (aX^/aX/vr) x (Xn/Xj,)* 
deviates from 1 by ~ 4%. We use this as an estimate in section [7] of systematic 
errors due to this effect. 

Finite lattice spacing 

Non-perturbative 0(a) improved clover fermions are employed in order to min- 
imise finite lattice spacing effects in the mass ratios determined here. Any effects 
are difficult to estimate if only one /3 (or a value) is available, so as a check 
we have performed some additional simulations at /3 = 5.80 (with an estimated 
a ~ 0.06 fm). The results are along the SU{3) flavour symmetric line (as three 
light quarks might show effects sooner than two light and one heavier quark). 
Again we have located (as described in [13]) the starting point on this line, mo for 
a path in the rris-mi plane leading to the physical point and then considered com- 
parable mass ratios for the nucleon, X^(m)/X^(mo) (against X^(m) / X^{mo)) , 
where we denote here a general point on the SU{3) symmetric line by m. In 
Fig. [9] we show these results. Both /3 values lie on a common line and show no 



1.5 



0.5 



P=5.50 
♦ P=5.80 

























0.0 0.5 1.0 _ 1.5 2.0 2.5 3.0 

x/( m)/x/(m„) 



Figure 9: X^(m)/X^(mo) against X^(m)/X^(mo) along the symmetric line. Square 
symbols are the /3 = 5.50 results and are given in [13) and [24j . while diamonds are 
the /? = 5.80 results [14]. (All results are on 32^ x 64 sized lattices.) The point where 
mo = m (ko = 0.12090 for /3 = 5.50 and kq = 0.12281 for /3 = 5.80) is denoted by a 
circle. 
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systematic lattice spacing dependence and so lie close to the continuum limit 
(certainly within the precision achievable here). 

B Tables 

We now give in Tables H] - 16] the baryon and pseudoscalar PQ mass results used 
in the analysis. Additional 32^ x 64 results along the line m = const, also used 
here are given in Tables 22 and 20 of [13] for the baryon and pseudoscalar meson 
particles, respectively. 









Kb 


M{aab) 





110000 





110000 


1.969(3) 





110000 





115000 


1.780(3) 





110000 





120000 


1.555(3) 





110000 





120512 


1.530(3) 





110000 





120900 


1.511(3) 





114000 





114000 


1.520(3) 





114000 





116000 


1.436(3) 





114000 





118000 


1.345(3) 





114000 





120000 


1.245(3) 





114000 





120900 


1.199(3) 





115000 





110000 


1.592(3) 





115000 





115000 


1.397(3) 





115000 





120000 


1.161(2) 





115000 





120512 


1.133(2) 





115000 





120900 


1.113(3) 





116000 





114000 


1.354(3) 





116000 





116000 


1.268(3) 





116000 





118000 


1.175(3) 





116000 





120000 


1.072(3) 





116000 





120900 


1.023(2) 





118000 





114000 


1.173(3) 





118000 





116000 


1.085(3) 





118000 





118000 


0.9887(25) 





118000 





120000 


0.8800(23) 





118000 





120900 


0.8267(23) 



Table 4: Results for the PQ baryon masses, M{aah) in lattice units at /3 = 5.50 from 
a 32^ X 64 lattice with sea quark masses at the symmetric point, i.e. ki = Kg = 0.12090 
and valence quarks Ka, Hh- 
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Kb 


M{aab) 





120000 





110000 


1.134(3) 





120000 





114000 


0.9726(26) 





120000 





115000 


0.9279(24) 





120000 





116000 


0.8812(24) 





120000 





118000 


0.7789(23) 





120000 





120000 


0.6588(23) 





120000 





120512 


0.6232(24) 





120000 





120900 


0.5945(26) 





120512 





110000 


1.081(3) 





120512 





115000 


0.8722(25) 





120512 





120000 


0.5945(25) 





120512 





120512 


0.5564(27) 





120512 





120900 


0.5247(30) 





120900 





110000 


1.041(4) 





120900 





114000 


0.8755(36) 





120900 





115000 


0.8295(32) 





120900 





116000 


0.7811(34) 





120900 





118000 


0.6739(33) 





120900 





120000 


0.5435(32) 





120900 





120512 


0.5031(35) 





120900 





120900 


0.4673(27) 



Table H] continued. 



1^1 


Kg 


M{sss) 


0.121040 


0.120620 


0.5265(16) 


0.121095 


0.120512 


0.5446(16) 


0.121145 


0.120413 


0.5682(13) 



Table 5: Additional results to Table 22 of [13] for the PQ baryon masses, M^^ = 
M{sss) in lattice units from a 32^ x 64 lattice along the line m = const. . 
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Kb 


M(ab) 





IIUUUU 





1 1 AAAA 

110000 


1 


2485 







110000 





1 1 r'AAA 

115000 


1 


0583 







110000 





1 OAAAA 

IzOOOO 





8351 


;4) 





1 1 AAAA 

110000 





1 AAAAA 

120900 





7909 


;io) 


n 

u 


1 1 AAAA 

IIUUUU 


A 




1 O A CT 1 O 

120512 





8100 


;6) 





"11/1 AAA 

114000 





11/1 AAA 

114000 





9436 


;3) 





11/1 AAA 

114000 





1 1 /^AAA 

116000 





8583 







11/1 AAA 

114000 





1 1 OAAA 

118000 





7664 


;3) 


U 


11/1 AAA 

114000 





1 OAAAA 

120000 





6669 


;4) 


u 


-i 1 A AAA 

114UUU 


A 




1 OAAAA 

120900 





6200 


7) 


(J 


lloUUU 





115000 





8593 


;3) 





1 1 AAA 

115000 


A 




1 OAAAA 

120000 





6202 


;4) 


u 


1 1 r AAA 

115000 


A 




1 OAAAA 

120900 





5720 


;6) 


U 


1 1 r AAA 

115000 


A 




1 on r 1 o 
120512 





5929 


^0 ) 


n 


1 1 fifinn 


n 


1 1 finnn 





7706 


;3) 





116000 





118000 





6754 


;3) 





116000 





120000 





5710 


;4) 





116000 





120900 





5213 







118000 





118000 





5752 


;3) 





118000 





120000 





4628 


[A) 





118000 





120900 





4077 







120000 





120000 





3342 


;4) 





120000 





120900 





2646 







120000 





120512 





2958 


;4) 





120512 





120900 





2174 







120512 





120512 





2534 


;4) 





120900 





120900 





1747 


;5) 



Table 6: Results for the PQ pseudoscalar masses, M{ab) in lattice units at /3 = 5.50 
from a 32'^ x 64 lattice with sea quark masses at the symmetric point, i.e. ki = Kg = 
0.12090 and valence quarks Kb- 
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